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and the dynamics is reduced by elimination of the electric field energy in the Hamiltonian. 
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I. INTRODUCTION 



Perturbation theory in the context of Hamiltonian dynamics has proved to be unquestion- 
ably useful in many contexts, ranging from celestial mechanics (e.g. to atomic physics 
(e.g. Q), to plasma physics (e.g. js]). The superconvergent expansions of Kolmogorov- 
Arnold-Moser theorem (e.g. [4]) and the techniques of adiabatic invariance (e.g. [5i]) all are 
aspects of perturbation theory in the Hamiltonian context. Although such techniques are 
well-developed and well-known for finite-dimensional systems, this is not the case for such 
perturbation theories for partial differential equations. This is particularly true for Hamil- 
tonian systems with noncanonical Poisson brackets of the form of those given in Uli for 
plasma systems. A main goal of the present paper is to provide tools for such perturbation 
theory using the Poisson bracket for Vlasov-Maxwell equations 0, l9l-[l2| in situations with 
a short time scale introduced by the presence of a strong magnetic field. 

Derivations of gyrokinetic theories have proceeded directly from the Vlasov-Maxwell equa- 
tions of motion as in the nonlinear development of {isl] , they have been based on Hamiltonian 
particle orbit perturbation theory that is lifted up to the kinetic level as in the linear de- 
velopment of [ij], or they have incorporated both particle orbit and kinetic perturbations 



to arrive at a nonlinear theory 



15|. (See 
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171] for review.) None of these procedures 



parallels that for finite-dimensional Hamiltonian systems that has historically achieved such 
great success. Consequently, none of these theories obtain an infinite-dimensional Hamilto- 
nian form as a consequence of their method of derivation. In fact, at present it is not known 
if nonlinear gyrokinetics has Hamiltonian form, the form possessed by all of the important 
systems of plasma physics when dissipative terms are neglected. (For review of Hamiltonian 
structure and techniques see [3, S Q)- 

To effect an infinite- dimensional Hamiltonian gyrokinetic-like perturbation theory re- 
quires a sequence of coordinate changes that involves both the dependent (field) variables 
and their arguments, which are independent variables from the point of view of the Hamil- 
tonian structure. This complicates matters significantly and care must be taken when per- 
forming transformations, most notably with the chain rule. Because the Vlasov-Maxwell 
theory has fields of mixed type, the electromagnetic fields depending on a space variable 
and the distribution function depending on a phase space variable, and because these fields 
are not a usual canonical set, the situation is further complicated. In the present paper 



the intricacies of this kind of transformation and associated chain rule are described, which 
enables the Hamiltonian perturbation theory. The techniques are then applied to obtain a 
simplified version of gyrokinetics (guiding center kinetics), which considers the presence of 
a conserved magnetic moment, as a first step for more general gyrokinetic reduction, e.g. by 



using the intrinsic coordinates developed in 19|, that will be considered elsewhere. 



The paper is organized as follows. In Sec. [TTl some preliminary material needed for the 
subsequent development is described. This is followed by four sections where several specific 
transformations are considered. It is shown how to lift these coordinate transformations, 
which are tailored to particle orbit dynamics, up to the level of fields by detailing how to 
transform the Vlasov-Maxwell Poisson bracket into the new coordinates. Lifting in this 



context is a natural relative of that treated in Ref. [12|], a companion paper that treats the 
lifting of microscopic particle dynamics up to the field level. Section IIIII considers magnetic 
coordinates, where the particle velocity coordinate is projected parallel and perpendicular to 
a space-dependent dynamic magnetic field. Because the new particle coordinates depend on 
the field, two chain rules must be considered: the usual function chain rule for phase space 
coordinates, where the field is assumed to be a given, and the chain rule for functionals 
which is needed for transforming the field theoretic Poisson bracket. Next, in Sec. llVi 
spherical velocity coordinates are considered. Here the velocity coordinates are chosen as 
the unit vector of the velocity (independent of the spatial coordinates) and a coordinate 
in one-to-one correspondence with the norm of the velocity. This transformation is a step 
closer to that needed for gyrokinetics where one introduces the magnetic moment, and it 
introduces a new feature in that the Jacobian determinant of the transformation is no longer 
unity. Thus, it provides a simple example for explaining how functional derivatives change 
when Jacobians change. In Sec. |Vl we turn to a more complete case where the change of 
coordinates depends only on the local value of the magnetic field, as a precursor to Sec. |VT] 
that considers the physically important situation, where the change of coordinates involves 
spatial derivatives of the magnetic field to arbitrary order, i.e., as given by Eq. ([T]) below. 
With the techniques of the previous four sections in hand, in Sec. I VI II we treat an example 
where the reduced coordinate is indeed the magnetic moment and explicitly transform the 
Hamiltonian form of the Vlasov-Maxwell equation into the new coordinates. Finally, in 
Sec. IVIIIt we conclude. 



II. PRELIMINARY MATERIAL 



From a general point of view, a main purpose of this paper is to transform the Vlasov- 
Maxwell Hamiltonian structure when the phase space variables (q, v) are changed to the 
following new coordinates that depend on the magnetic field and all of its derivatives: 

q = q, v = v(q,v;B,VB,...) (1) 

For the noncanonical Hamiltonian structure of Vlasov- Maxwell dynamics, the observables 
are the set of all functionals of the magnetic field B(q), the electric field E(q), and the phase 
space density /(q, v), where the time variable has been suppressed. The Poisson bracket is 



imii\: 



{F,G} = jd\d\f \Ff,Gf] (2) 

+ e jd\d^v f (Ge • d^Ff - Fe • d^Gj) 

+ Jd^q (Fe ■ V X Gb - Ge ■ V X Fb) , 

where subscripts are used for functional derivatives, Ff := 6F/6f, Fe := 5F/5E, etc., and 
the particle bracket is [/, g] = V/ • d^g — Vg ■ d^f + eB • d^f x d^g, with V/ = df/dq^ and 
d^g = dg/dw. For the sake of simplicity physical constants have been scaled away as usual, 
but a dimensionless charge variable e that indicates the coupling term has been retained 



l^see for a dimensional form of this bracket). The variable e becomes the charge ratios 
when ([2]) is generalized by summing over multiple species. 
The Hamiltonian functional is 

i/[E,B,/] = l Jd'qd\\\vrf + l Jd'q iWEf + iBf) , (3) 

which is the sum of the kinetic energy of the plasma and the energy of the electromagnetic 
field. The relativistic model is obtained by replacing ||v|p in the kinetic energy term with 
a/1 + ||vP, where in the latter case v is the scaled relativistic momentum. The coupling 
between the plasma and electromagnetic field is included in the noncanonical Poisson bracket 
([2]). The Hamiltonian ([3]) together with the Poisson bracket generates the motion through 
Hamilton's equations expressed as 

F = {F,H}, 
4 



for any observable F. In particular, if F denotes the field variables the bracket induces 
Maxwell- Vlasov equations as follows: 

B = {B,i/} = -V X E, 

E = {E,/7} = V X B -e jdh /v, 

/ = {/, /f} = -V V/ - e (E + V X B) ■ 9v/ . 

As noted in Sec. HI in order to transform the Hamiltonian structure to facilitate the sepa- 
ration or removal of fast time scales (as in oscillating-center, guiding-center, and gyrokinetic 
theories) care must be taken because such a change of coordinates involves both the depen- 
dent and independent variables, i.e., the spatial observation points of the field. A simple 
case of this is treated in the next section. 



III. LIFTING WITH MAGNETIC COORDINATES 

As a first case of lifting, consider velocity coordinates based on a decomposition of the 
velocity using the magnetic field. This transformation of the spatial coordinate is unchanged, 
but the velocity v is transformed as follows: 

V = v(v; B) = v(t;||, v_l; B) = b^n + v_l , 

where b = B/||B|| is the unit vector the direction of the magnetic field, 

= b • V 

is the (scalar) component of the velocity parallel to the magnetic field, and 

= V — bb ■ V = ■ V 

is the (vectorial) component of the velocity perpendicular to the magnetic field, with 

h:=I-hh (4) 

being the orthogonal projector onto the plane perpendicular to the magnetic field. 

There are two chain rules to consider: that for functions, considered next, and that for 
functionals, such as the energy expression of ([3]), which will follow. 
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A. Function chain rule 



The transformation of the field Poisson bracket of ([2]) requires the transformation of the 
particle bracket, 

|/.9l = ^-|^-^# + eB.f|ix|»V (6) 



(9q (9v 9q dv \dv dv 

into the new coordinates, (q, v) — )■ (q, f||,v^). The following abbreviations are convenient: 

V := — , a = — , d\\ = —— , di 



dc{ ' dqi ' dv\\ ' d\± ' 

Note the last operator acts only in the plane perpendicular to B, which implies the following 
properties: 

d±f-h = dj and h.dj = 0. 
Total variations of /(q, v) = /(q, t>||, vj_) are given by 

5/ = 1^ ■ 5q + |{ ■ 5v = V/ ■ (5q + + dj ■ , (6) 

while variations of the initial and final coordinates are related by 

5^11 = b ■ (5v + {5q ■ Vb) ■ V , 
5v_L = /j_ ■ (5v — 6l± ■ V 

= Jx ■ (5v - (5q ■ Vb) (b ■ v) 

-b(5q-Vb)-v. (7) 

For the function chain rule the field B is assumed to be a fixed function with the coordinates 
(q, v) changing. 

Inserting ([7]) into ([6]) implies the chain rule relations 

1^ = b dj+dj-h = b + dJ, (8) 

|Z = dJ + (v ■ d,h) - (b ■ v) dJ ■ d,h , (9) 

and using ([8]) and (Ej) in ([5]) gives the particle bracket in the magnetic coordinates 

[/,^] = b-(V/9||^-V^a||/) (10) 

+ a- {d±g d\\f - d±f d\\g) 

+ dj-h- d^g + eB ■ {dJ x d±g) , 
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with 

tti — V ■ dih + (b • v) b • dhi and = (b • v) (^dihj — djhi^ . 

In all these relations, recall that d±f = d±f ■ I±. This is important because, for instance, 
the component of V^, a or b parallel to b are non-zero, but vanish when contracted with 
dj. 



B. Jacobian 

In general care must be taken with the Jacobian determinant when defining functional 

derivatives, but here the Jacobian is unity 

<9(q,v) 

This follows because rotations have unit Jacobians and at any time there exists a rotation 
to a cartesian coordinate system with one of the v axes aligned with b. Thus 

dz :— (fqd^v — d^qdv\\d'^v± =: dqdv . 

Because the volume integral is ultimately independent of how it is calculated, dz can be 
assumed to be independent of B, e.g. when calculating functional derivatives with respect 
to B, the topic considered next. 



C. Functional chain rule 

For the functional chain rule, the transformation of the fields must be made definite. 
Here, 

E(q) = E(q) , B(q) = B(q) , 
/(q, v) = /(q, ^11, v±) = /(q, b • V, • v) 
= /(q, bT;|| + v±) , 

where now the coordinates (q, v) are fixed and the field b varies. 

Variation of a transformed functional, F[f, B, E] = F[f, B, E], gives 

SF = Jdz FfSf + Jdq (Fb -SB + Fe- SE) 

^ Jdz Ff-Sf + jdq (Fg -SB + Fe- SE) . (11) 



With the variations of the initial and final fields related by 



5E = 5E, SB = SB, and 5f = Sf + diifi^vSh) + d±f-SI±-v 



[12] 



expressions relating functional derivatives of new and old variables can be obtained. Using 

1 



B 



SB-dJ^ 



and after some work the last equation of (1121) becomes 

Inserting this and the other two equations of f lT2|) into ffTTl) . and then equating coefficients, 

gives the functional chain rule relations 

6F _6F 5F _ SF 
If ~~6f' 5E ~ 5E ' 

(13) 



6B = 6B^mJ "W^-^^^ 



where 



d;:=^f±d\\-v\\d^. (14) 
Finally, the Maxwell- Vlasov bracket expressed in these magnetic coordinates is 



{F,G} = jdzf[Fj,Gf] 

+ e Jdz /{Ge- d^Ff - Fe ■ d^Gj 



(15) 



+ d^q i^E-Vx 



Ge + 



IBI 



dv Gf dlf 



-Ge ■ V X 



IBI 



dv Ff d*J 



where the 'bars' have been dropped, [, ] means the bracket of (|5]) rewritten in the new 
coordinates as (fTOj) . and = h d\\ + 9^ is a shorthand as in ([8]). Note, 9*f ^ = 0. 



IV. LIFTING SPHERICAL VELOCITY COORDINATES v = y v 



Now turn to the new coordinates considered for intrinsic gyrokinetics (used in 19] ) , which 
changes only one of the velocity coordinates to get the magnetic moment. The two other 
velocity coordinates are usually chosen as the unit vector of the velocity. So, a preliminary 
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change of coordinates consists in adopting spherical coordinates for the velocity space: v = 
Vir where V := ||v|| G M+ is the norm of the velocity and v := v/||v|| G is the unit 
vector of the velocity. This transformation is considered in this section, but later the change 
y — 7- yU will be considered. 

The transformation v ^ (v, V") is clearly invertible. For the chain rule the following are 
needed: 

= 5v 

5V = V ■ 5v and 5v = Xj_ ■ — , 

where 

X± = I — vv 

is the orthogonal projector onto the plane perpendicular to the velocity. Note Xj_ is different 
from the magnetic projector I± of (jl]) used in Sec. IIIII 

As in Sec. IIIII the above are used to calculate the function chain rule, giving 

df Idf = df ^ 

Inserting f|T6|) and f|T7j) into ([5]) and, after some manipulations, the particle bracket ex- 
pressed in spherical coordinates is obtained 

[/, g] = ^ (V/ ■ ■ d^g -Vg-X^- d^f) 
+ V ■ (V/ dyg - Vg dvf) 

+ . (^Ovf d^g - d^f dyg) , (18) 

where, for convenience, the 'bars' have been dropped and the abbreviations 

%=--d.f and %=--dvf. 

have been employed. 

Turning to the functional chain rule, notice that the change of coordinates does not 
depend on the fields, but the Jacobian for this special case is not unity 



dz = V^dVdQd^q = JdVdnd^q =: Jdr]d^q =: Jdw , 
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(19) 



because the integration measures are changed from d'^v and dz to drj and dw, which are 
defined by relation f[T9|) . 
Thus, as above, 

5F = jdz FfSf + j d^q (Fb -SB + Fe- 5E) (20) 
= jdwFf6f + jd\ (Fb ■ + Fe ■ . 
Inserting f|T9|) into fl20|) gives 

Ff = J-'Ff-, Fb = Fb, and F^ = F^ . (21) 

Note, in (12T!) the new functional derivative is defined with respect to the bare measure dw. 
So, the first term of the Maxwell- Vlasov bracket transforms as 

{F, G}i := Jdzf [Ff, Gf] = jdw Jf [j-'Fj, J-^Gj] = {F, G}, , (22) 

with the bracket of the second equality above given by f|T8l) . 

The basic identity for this bracket with Jacobians, which replaces the usual ^f-g-h^ iden- 
tity for canonical brackets Jdz f[g,h] = — J dz g[f, h], is the following: 

jdw Jf[j-'g,h] = - jdw g [/, h] . (23) 

In terms of the bare measure 



i/(u,') 

The bracket of (1221) with (!23l) and (12^ produces the correct equations of motion for the 
Vlasov-Poisson system. 

Now consider the coupling terms of the bracket 

{F, G}2 ■.= ejdzf (Ge ■ d^Ff - Fe ■ d^Gf) 

= e jdw Jf (Ge ■ d^J-^Fj - Fe ■ d^J-^Gj) , 

where is a shorthand for the expression of f|T6l) . When generating Maxwell's equations, 
the Hamiltonian gives 

which gives the correct expression for the current density 3 = j dr] jT/v. 
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Finally, the pure field terms of the Maxwell- Vlasov bracket are unchanged and, thus, the 
Maxwell- Vlasov bracket in these spherical coordinates becomes 



where the 'bars' have been dropped, and [ , ] means the bracket of rewritten in the new 
coordinates as (ITS]) . 

V. LIFTING WITH LOCAL DEPENDENCE ON B 

To include the magnetic moment in the coordinates, the next step is to investigate the 
coordinate transformation V i-^ A, where A is a coordinate in one-to-one correspondence 
with the coordinate V of Sec. HVl but in this section it is assumed to have local dependence on 
the magnetic field, i.e, it depends on B but not its derivatives. Explicitly, the tranformation 
is (q, V, v) O (q. A, w) where 



Clearly, invertibility requires Va '■= dV/dA ^ 0. Since the first two equations above are 
identities, eventually v will be used for w and q for q. 
The Jacobian for this transformation is now 



which define the Jacobian J' and the integration measures drj and dw. Note that these are 
not the same as those of Sec. llVt even though the same symbols are used. Furthermore, 
now depends on B and, hence, q. Also, dQ contains a portion of the Jacobian from cartesian 
coordinates, but one that is independent of q. 

Now the chain rule is effected on functions analogous to ([8])- ([9]) and (fT6|) - (fT7|) and on func- 
tionals analogous to (|T3|) and ( pTj) . Varying /(q, V, v) = / (q. A, w) in the label (coordinates) 




q 



q, v = w, and V = V{A,w,'B). 



dz = V^dVdQd^q = V^VAdAdQd^q 
= JdAdVld^q =: Jdrjd^q =: Jdw , 
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dependence, and then equating as above, gives 

dl^dl_VB^m. dl 

(9q (9q Va <9q dA ' 

dl^J_dl 

dV Va dA ' 

9/ _ 9/ _ J_dV_ dl_ 

dv dw Va dv dA 

Inserting ( l26ll and (|27l) into (fT6l) gives the chain rule on functions 



while ( l25l) gives 



*^ ~ dv~ V \dw ~ VaOA dw) ' VaOA 



(9q 9q Va 9q Va 



Then, inserting ( l28l) and ( l29l) into ( |T8l) gives the following complicated expression for 
particle bracket [ , ] in the new coordinates: 

[/, g] = V./ ■ D,-g - V.g ■ Dj + eB ■ [dJ x D^g) 
= ^ (V/ ■ ■ d^g -Vg-t^- d^f) 

+ ^-{Vf dA-g-VgdAf) 

+ ^%^-{V-gdAf-VfdA-g) 

+ ^VBi-l^- {d^f dA-g - d^-g dAf) 

+ ^ ■ (d^f ■ X (^d^g ■ Xx) 

+ ^ ^ ■ [d^g dAf - 9w/ dAg) 
V Va 

X fx, -d^V 



V^Va 

■ Xj_ ■ {d^g dAf - 9w/ dAg) ■ 



Now consider the functional chain rule as above, 



5F = jdzFf5f + jd^q {Fb ■ 5B + Fe ■ 5E) 
= jdwFf5f + jd\ (Fg ■ + Fe ■ 
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Functionally varying /(q, V, v) = /(q, A, w) gives 



(32) 



while 5B = 5B and 5E = (5E. Whence, upon substitution of fl32|) into fl3T|) . the chain rule 
on functionals is obtained, 

5F _ 1 5F 
Jf~jJf' 
SF _ 5F 

6E~6E' 
SF _ SF 

sb~sb' 



dA df SF 



(33) 



dB dA Sf ' 

where the last expression of ( l33ll can be written in a more convenient way as 

S£^SF r \^ dls_F 

SB SB J Va dA Sf ■ 

This follows from 

dA _ 

dB ~ ~l4 ' 

which comes about because the change in A induced by a change in B at fixed V and w, 
satisfies = SV = VaSA + Vs^SBi. 

Finally, the Maxwell- Vlasov bracket in the coordinates (q. A, v) is given by 



{F,G} = jdvd^q Jf [j-'Ff,J-'Gf] 

+ e jd7]d\ Jf (Ge ■ D,J-'Ff - Fe ■ D,J-'Gf) 



(34) 



+ d^q Fe- V X 



dr] 



Vb dfSG 



Fb+ dr] 



Va dA Sf 

Vb df sf 



- ' Va dASf 

where the particle bracket [ , ] is given by ( 15U]) . is the operator defined by 
bars have been dropped. 



and the 



VI. LIFTING WITH NONLOCAL DEPENDENCE ON B 

In order to include the physical coordinates where A is the magnetic moment fi, the last 
step is to consider the case where the coordinate transformation involves derivatives of the 
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magnetic field. This is important because perturbative reductions, such as those based on 



Lie-transforms 



3-21 



2l| or mixed variable generating functions 22|], often involve derivatives 
to arbitrary high order in the fields. 

So, a more general transformation to new coordinates (q, V, v) ^ (q. A, w) is considered: 

q=q, v = w, and V = V[A,w,B], (35) 

where now l^[yl,w, B] means a transformation that depends on B and, possibly, all its 
derivatives. Clearly, invertibility requires Va '■= dV/dA ^ 0, as in Sec. El Since the first two 
equations above are identities, as before eventually v will be used for w and q for q. 
The Jacobian for this transformation is again 

dz = V'^VAdAdVtd^q = JdAdnd^q =: Jdrjd^q =: J dw . 

but now J depends on q through B and its derivatives. 

For the chain rule on functions or functionals, /(q, V,v) = /(q, A, w) is varied as in 
Sec. lYl and all terms are the same as before, except some slight changes in the relations 
involving derivatives with respect to the magnetic field. Indeed, the Frechet derivative with 
respect to B is now a differential operator, and care must be taken with the order of terms. 
For instance, relation fl25p becomes 

dq dq OA Va ' dq ' ^ ' 

where Vb. is now a differential operator acting on dBi/dq. Formulae (I29l)-( l30|) must be 
changed accordingly. 

As for relations ( l32l) - (!33|) . variation is performed slightly differently this time as follows: 

Sf = 6f + fv V^-6B, 

where is the Frechet derivative operating on 5B. Thus the chain rule for functional 
derivatives gives 



5B - 5B + y "^"^ [dV 6f) " 5B + y "^"^ Ua OA 



where the adjoint f is done with respect to dw. 
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Finally, the Maxwell- Vlasov bracket in these coordinates becomes 
{F,G} = Jdr^d'q Jf [j-'Fj,J-'Gj] 
+ e jdrid^q Jf (Ge ■ D,J-'Ff - Fe ■ D,J-'Gf) 
(F,.Vx[Ge + /*,r^(||^)' 
-Ge ■ V X 



(37) 



VII. AN EXAMPLE USING THE MAGNETIC MOMENT 



With the transformed bracket (l37Il , the first thing to be checked is whether the dynamics 
agrees with the conservation of the magnetic moment, when appropriate, since this is what 
suggested the reduction in the first place. To this end, suppose the coordinate A is the 



magnetic moment, A := /i(q, v), which to lowest order is given by A 



|B||. To get 



a true conserved quantity, small corrections must be added to all orders in the Larmor radius, 
including derivatives of all orders in the magnetic field [l6|,ll9|. Thus, A = || vxp/||B|| +0(e) 
is defined as solution of the following equation 

= /i = V ■ V/i + ev X B ■ . 

At the field level, the conservation of the magnetic moment corresponds to the conserva- 
tion of the functional 

M := jdz ffi 

for any particle distribution /. In the transformed coordinates, this is 

M := Jdw Jf^i. 

To investigate the conservation of M, note that a static magnetic field corresponds to 
elimination of the electric field term in the Hamiltonian functional, since this eliminates the 
V X E term in the Maxwell-Faraday equation. In this case 



M = {M, H} = jdr]d'q Jf [fi, J'^Hj] 

= i JdridS Jf (V,/i ■ D,V^ + eB ■ D^fx x D,V^) 
d\d^q / (v ■ V/i + V X eB ■ d^^i) = , 
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as was expected. 

Accordingly, the transformed bracket fl37|) is expressed in coordinates adapted to the 
conserved magnetic moment. As is ususal in gyrokinetics, the electromagnetic field dynamics 
spoils the conservation of the magnetic moment. This is why the feed-back of the plasma 
dynamics onto the electromagnetic field dynamics needs to be restored as a perturbation, 
i.e., a perturbed magnetic moment must be defined that is conserved 

Consider now the transformed Maxwell- Vlasov equations of motion generated by the 
bracket (!37|l . In this bracket, most of the terms are actually identical to those of the initial 
bracket ([2]), even though their formal expressions look different because they are expressed 
in the reduced coordinates (q, A, w), e.g. through formulae (!28|) and (1301) . The only new 
terms are 



and one obtained by permuting F and G (and with a minus sign for bracket antisymmetry) . 

In the equations of motion, this new bracket term generates an additional term in 
Maxwell-Ampere equation, viz. 



At first glance this additional term looks like a new magnetization current. But, one must 
remember that the usual V x B term has itself another additional contribution V x SHf^in/ SB, 
because in the reduced variables, the plasma kinetic energy depends on the magnetic field 
Hkin '■= JdwJ'fV'^/2 that is not constant in B (both because of J' and V). And, it turns 





(38) 
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out that this last additional contribution exactly cancels the "magnetization" term in (138|) : 



6BM 



I Jdw f{dAV-V'')^5{q-x.) 
JdwdAf{V')^6{ci-^) 
JdwdAf{V')^V^S{q-^) 

-\ jdw OaI JVVa V^b{<q - x) 

- ^dw Ba] § V^K^. - x) 

-JdWi^^dAf) . 



This cancellation was to be expected, since the electric field E = E is not affected by the 
change of velocity coordinates, and the current term has not been changed either, but only 
expressed in the new variables: 

-e JdT] Jf {J-'Hf) = -| Jd^'v f 9v {J-'JM^) = -J . 

Finally, the additional term in the transformed bracket (I57|) generates another additional 
term in the equation of motion: the dynamics of the Vlasov phase space density / has an 
additional force term 

1 df 

This term is not cancelled by any other term. It can be rewritten as 

which is exactly the expected contribution when applying the chain rule for the time deriva- 
tive of the transformed fields. It comes about because the change of coordinates is time- 
dependent when the magnetic field is not static. 



VIII. CONCLUSION 

In summary, in this paper techniques for transforming the Vlasov-Maxwell Poisson 
bracket to new coordinates, when the transformation law mixes dependent and independent 
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variables, have been developed. Four transformations were considered, each of which con- 
sidered a new feature needed for understanding the more general transformation of fl35l) . In 
Sec. mil a transformation that mixed the independent velocity variable with the magnetic 
field was considered and the associated function and functional chain rules were described. 
In Sec. \IV\ spherical velocity coordinates were treated and here it was seen how a nontrivial 
Jacobian determinant influences a transformation. In Sec. |V] a class of transformations 
that mixes the dependent and independent variables by having dependence on B and in 
addition possesses a nontrivial Jacobian was considered. Finally, in Sec. |VT1 the nonlocal 
transformation of was effected, the most general transformation of this paper that 
results in the transformed noncanonical Poisson bracket of (P7|) . This final form of the 
Poisson bracket was seen to contain additional terms that appear to be magnetization-like 
contributions. However, these bracket terms were shown to produce no magnetization term 
in the equations of motion, since the electromagnetic fields are not affected by the change 
of field coordinates. Only the dynamics of the Vlasov density obtained an additional term, 
a term that results from the change of field coordinates being time-dependent through B. 

The transformations of Sees. IIIIHVII paved the way for the simple example of Sec. IVIII 
Here the dynamics was reduced by dropping the electric field energy from the Hamiltonian, 
resulting in the magnetic moment being conserved by a reduced dynamics that must have a 
static magnetic field. However, when restoring the feed-back of the plasma dynamics onto 
the electromagnetic field dynamics, the magnetic moment was seen to be no longer conserved 
and must be perturbatively changed to be conserved. 

In all the cases considered, the lifting was eased because the change of coordinates only 
concerned a new particle velocity that depends on the magnetic field, but no change was 
made in the spatial coordinate. If Eq. ([T]) is generalized by adding dependence on the electric 
field and all its derivatives, then results similar to those presented are immediate. However, 
if the new spatial variable has velocity and field dependence, then the situation becomes 
considerably more complex. Such transformations are of interest for some oscillation-center, 
guiding-center, and gyrokinetic theory development, and the same methods of function and 
functional chain rule can be used, but some additional effects will show up, e.g., non-zero 
polarization and magnetization terms like those of 12|]. Details of the magnetic moment 
reduction will be given in 19| and more general lifting will be considered in a future publi- 
cation. 
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